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ABSTRACT 

We estimate the observed distribution of chirp masses of compact object 
binaries for the gravitational wave detectors. The stellar binary evolution is 
modeled using the StarTrack population synthesis code. The distribution of the 
predicted "observed" chirp masses vary with variation of different parameters 
describing stellar binary evolution. We estimate the sensitivity of the observed 
distribution to variation of these parameters and show which of the parameters 
can be constrained after observing 20, 100, and 500 compact object mergers. 
As a general feature of all our models we find that the population of observed 
binaries is dominated by the double black hole mergers. 

Subject headings: binaries: close — gravitational waves 

1. INTRODUCTION 

Compact object mergers are one of the most promising sources of gravitational waves for 
the ground based interferometric detectors like LIGO (Abramovici et al 1992) and VIRGO 
(Bradaschia etal 1990). So far most of the theoretical papers on the properties of these 
sources related to the gravitational wave detections have been concentrated on calculation 
of the predicted rates (Narayan etal 1991, Phinney 1991, Kalogera etal. 2001). In this paper 
we wish to address another aspect of the gravitational wave detection, i.e. the distribution 
of observed masses of the compact objects. 

Stellar mass compact object binaries shall be detected during the inspiral phase, while 
the consecutive merger and ringdown phases will most likely have lower signal to noise ratios. 
During the inspiral phase the motion of the binary components and also the wave form is 
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governed by the chirp mass M = (m 1 +m 2 ) 2/5 ( y m 1 m 2 ) 3/5 (Peters and Matthews, 1963). The 
waveform will depend on the individual masses of the binary components m\ and m 2 when 
the post Newtonian effects are taken into account. However, the analysis of the inspiral 
phase alone shall not suffice to determine if a binary contained a neutron star or a black hole 
without the prior knowledge of the neutron star maximum mass. A careful modeling of the 
signal may yield the individual masses of the objects, however, the chirp mass will be the 
primary observable for the compact object mergers (Cutler & Flanagan 1994). 

We use the StarTrack population synthesis code (Belczynski, Kalogera & Bulik 2002) 
to calculate the distributions of compact object binary masses, and we present these calcu- 
lations in § 2. In § 3 we estimate the number of merger detections required to distinguish 
between different models of stellar binary evolution. Finally, § 4 contains the conclusions 
and discussion. 



2. Distribution of the chirp masses 

The StarTrack binary population synthesis code is described in detail in Belczynski, 
Kalogera & Bulik (2002). One of the important features of the code is the possibility to 
conduct parameter study of a given property of the population of binaries, i.e. to estimate 
the dependence of the result on each of the parameters used to describe the stellar and 
binary evolution. The models used in this paper are listed in Table 1. We first use the 
standard model A results to present the intrinsic distribution of the chirp masses. This is 
shown in Figure 1. The distribution shows a clear peak at low chirp masses 1M Q M. < 2M Q 
which is due to the double neutron star systems. The mixed (BH-NS) systems populate 
the intermediate region, while the chirp masses of the BH-BH binaries extend up to above 
1OM . 

In order to estimate the observed distribution of the chirp masses of compact objects 
one has to take into account the sensitivity of the gravitational wave detectors to signals 
from mergers of different binaries. The calculations of the signal to noise ratio (Finn and 
Chernoff 1993, Bonazzola and Marck 1994, Flanagan and Hughes 1997) show that the sam- 
pling distance in the first approximation is a function of the chirp mass only: D oc A4 5 ^ 6 . 
The additional corrections due to limited sensitivity window of the detectors have been cal- 
culated by Flanagan and Hughes (1997) and amount to approximately 10% for the binaries 
with the total mass below 18M for the initial LIGO, and less for the advanced LIGO. In 
this paper we neglect these corrections. The distribution of the expected observed chirp 
masses can be calculated using Monte Carlo method. We assume that the Universe is uni- 
formly filled with merging binaries, and for each merger we estimate the signal to noise ratio 
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Table 1: Description of different population synthesis models for which the distributions of 
mass ratios have been found. 

Model Description 

A standard model described in Belczynski, Kalogera, Bulik 

(2002), but with T Hubble = 15Gyrs 

Bl,7,ll zero kicks, single Maxwellian with a = 50, 500, kms -1 , 

B13 Paczynski (1990) kick with V k = 600km s" 1 

C no hyper-critical accretion onto NS/BH in CEs 

El-3 a CE x A = 0.1,0.5,2 

Fl-2 mass fraction accreted: f a = 0.1, 1 

Gl-2 wind changed by / w i n d = 0.5, 2 

J primary mass: oc Mf 2 ' 35 

Ll-2 angular momentum of material lost in MT: j = 0.5, 2.0 

Ml-2 initial mass ratio distribution: $((/) oc q~ 2 ' 7 ,q 3 

O partial fall back for 5.0 < M co < 14.0 M 

S all systems formed in circular orbits 

Zl-2 metallicity: Z = 0.01, and Z = 0.0001 
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Fig. 1. — The intrinsic (galactic) distribution of the chirp masses in the framework of model 
A (left panel), and the distribution of the expected observations (right panel). The solid line 
corresponds to the NS-NS mergers, the short dashed line represents the NS-BH mergers and 
the dashed line stands for the BH-BH margers. The sum of the three distribution in each 
panel is normalized to unity. 

in the detector. We model the population of merging binaries assuming a continuous star 
formation rate. The result is shown in the right panel of Figure 1. One can note that these 
distribution could also be obtained analytically by multiplying the distributions of Figure 1 
by the volume oc Ai 5 ^ 2 and normalizing it. In this plot the BH-BH systems are now the 
dominant contribution of the distribution. This is due to the fact that the sampling volume 
for the BH-BH binaries more than 100 times larger that that for the NS-NS systems, which 
easily compensates for the lower merger rate of the BH-BH binaries. 



3. Expected observations 

Let us now address the following questions: are the distributions of observed chirp 
masses expected in the framework of alternative models different? If so, are these differences 
significant? We simulate the distributions of chirp masses in the expected observations with 
the binary populations obtained from the set of models of Table 1 similarily as we have 
done above for the model A. We present the results in Figure 2. Different stellar evolution 
models lead to drastically different distributions of the chirp masses in the expected obser- 
vations. Various parameters describing stellar evolution affect the distribution of observed 
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chirp masses in several ways. Changing the kick velocity distribution (models B) alters the 
ratio between the number of the neutron star binaries and the back hole binaries. Other 
models change the maximal masses of the black holes produced. This is especially clear in the 
case of models G where the stellar winds are varied by a factor of two upwards (model G2) 
and downwards (model Gl) We note that the shapes of these distributions do not depend 
on the sensitivity of a detector. 

In order to verify if the differences between the distributions are significant we turn to 
a simulation of a finite number of merger observations. We assume that the true stellar 
evolution goes through one of the models of Table 1. We then simulate the observations of 
a given number of mergers (we use 20, 100 and 500 mergers) and for each such simulated 
observation we verify using the Kolmogorov Smirnov (KS) test if we can reject a hypothesis 
that the stellar evolution is described by model A. This allows to test the sensitivity of 
the shape of the distribution of expected chirp mass observations to the underlying model 
parameters describing stellar evolution. For each number of merger observations we repeat 
this test 10000 times to obtain a distribution of KS-test probabilities and find the lowest 
probability that appeared in the top one percentile of this distribution. We can now set a 
detection confidence level, say at 10~ 5 and compare each probability with this value: if it is 
higher we conclude that this particular model cannot be distinguished from model A with a 
given number of merger observations, while a smaller number means that this model can be 
distinguished, and that some constriants can be iposed on the particular parameter through 
which this model differes from model A. We present the results of the test in Figure 3. 

Figure 3 presents a measure of sensitivity of the expected observed distribution of chirp 
masses to the parameters describing stellar evolution. One can see from Figure 3 that even 
observations of a small number of mergers (open circles correspond to 20 mergers) yield 
highly significant results for models El, G2 and O. The reason for that can be is clear from 
Figure 2. These are the models for which the maximal chirp mass in the population is 
significantly lower than that for model A. Model G2 population (i.e. with increased stellar 
winds) contains hardly any black holes. In general we see that these observations are very 
sensitive to the value of maximum mass of stellar black holes in the population. Model Gl 
(with decreased stellar winds) which allows for formation of black holes binaries with chirp 
masses up to 16 M will stand out with less than a hundred merger observations. 

With a larger number of merger observation (stars in Figure 3) correspond to 100 merger 
detections) more parameters can be constrained. Some constraints can be obtained for the 
value of common envelope efficiency ac^A (models E: model E2 is similar to A). Other 
parameters describing mass transfer events like the mass fraction accreted (models F), and 
the amount of angular momentum lost (models L) shall also be constrained. Moreover also 
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the metallicity of the progenitor stars may influence the observed distribution at a significant 
level (models Z). 

Constraining the initial mass ratio distribution (models M) will require an even higher 
number of merger detections: only for the case of 500 observed mergers the differences 
become significant. Models C (no hypercritival accretion onto a compact object), J (initial 
mass function slope), and S (systems circular initially) lead to very small differences in the 
observed distribution of chirp masses. 

Models B where the kick velocity distribution is varies begin to show significant dif- 
ferences only with a large number of observations. Changing the kick velocity distribution 
affects strongly the absolute rates (Lipunov Postnov Prokhorov 1997, Belczynski Bulik 1999), 
and the ratio of double neutron star mergers to the double black hole mergers (Belczynski, 
Kalogera and Bulik 2002). 

4. CONCLUSIONS 

We have applied the stellar population synthesis models to simulate the distribution 
of observed chirp masses in the gravitational wave detection of stellar mergers. We find 
that the population of observed mergers is dominated by the the black hole - black hole 
binary mergers. In most models double black hole mergers constitute more than 90% of 
the observed events. The exception is model G2, in which the formation of black holes is 
suppressed because of increased stellar winds. The shapes of observed distributions of chirp 
masses vary considerably for different models of stellar binary evolution. 

We simulate the observed distributions of chirp masses in the framework of various stellar 
evolution models and estimate the sensitivity with which these parameters can be estimated 
from a given sample of observed mergers. We find that there is large number of parameters 
that can be constrained given a set of measured chirp masses. The main and immediate 
constraints come from the fact the the observed population seems to be dominated by the 
highest mass black hole binaries. Thus even a small set of observations yields constraints on 
the maximal mass of merging black hole binaries. A larger set of observations will lead to 
constraints on the evolution of high mass binaries. 

In our simulation we use a simple statistical tool: the Kolmogorov Smirnov test. Given 
a set of real observations with some measurements of individual masses of coalescing stars, 
one could use a more sensitive tool like the maximum likelihood method. However, even 
with such simple statistic as used here we can show the general properties of the expected 
observations, and demonstrate the sensitivity of the observed distributions to different model 
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parameters. 

we note that consideration of the distribution of observed masses will lead to stricter 
constraints than consideration of just the observed rates. The theoretical calculation of 
rates involves estimating a number selection effects and calibrating with other sources which 
leads to several uncertainties. The calculation of the observed distribution of chirp masses 
is free from such uncertainties because a distributions essentially equivalent to considering 
the ratios of the number of mergers of different type, and all the normaliztion factors that 
enter the rate estimates do cancel out. 

Finally, we have to mention several effects not taken into account in this paper. A more 
detailed calculation must include the consideration of the effects of lifetimes of binaries of 
different type. Belczynski, Kalogera and Bulik (2002) have shown that the typical lifetimes 
of double neutron star binaries are much smaller than the black hole binaries. The effects 
due to changing of the star formation rate with the redshift will affect the observed popu- 
lation of merging binaries. When considering the advanced detectors sensitive to mergers 
at cosmological distances one also needs to take into account cosmological effects: the fact 
that the true measure quantity is the redshifted mass (1 + z)M., and also the change of the 
observed volume with the redshift (here we have assumed the Euclidean geometry). These 
issues will be considered in a separate paper. 

We thank Vicky Kalogera for comments on this project. TB thanks the hospitality of 
the Theoretical Astrophysics Group at Northwestern University. We acknowledge support 
from KBN through grant 5P03D01120. 
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Fig. 2. — Distributions of the expected observed chirp masses in the framework of models 
listed in Table 1. For clarity each distribution is shifted up by a factor of ten. 
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Fig. 3. — The significance of rejecting model A. Open circles correspond to observations of 
just 20 mergers, stars to 100 mergers, and filled circles to observations of 500 mergers. The 
symbols with an arrow denote the case when significance is off the scale. 



